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\-{ • Abstract 



In this work, we employ the effective vertices for interaction between 
diquarks (scalar or axial-vector) and gluon where the form factors are 
derived in terms of the B-S equation, to obtain the potential for baryons 
including a light quark and a heavy diquark. The concerned phenomeno- 
logical parameters are obtained by fitting data of S^*-*— mesons instead 
of the heavy quarkonia. The operator ordering problem in quantum 
mechanics is discussed. Our numerical results indicate that the mass 
splitting between B^/2(^)i Bi/2{V) and Bi/2iS) is very small and it is 
consistent with the heavy quark effective theory (HQET). 



I. Introduction 



It is well known that the heavy flavor physics can be different from the world where only 
light flavors are involved in many aspects. Since the heavy flavor can serve as a static color 
source as Mq » ^qcd-, an extra symmetry SUf{2) ® SUs{2) exists |j^] which can attribute 
all non-perturbative QCD effects into one or a few form factors and make the hadronic matrix 
elements evaluation much simplified. On the other hand, the diquark structure in baryons causes 
interests of many theorists of high energy physics Q. The possible diquark structure in nucleons 
has been studied in the non-relativistic QCD-based potential quark model Q. Even though as 
pointed out, a spin-0 diquark structure may exist in nucleons, one has reason to doubt its 
validity. Because the two light quark which are supposed to constitute a diquark are relativistic 
and dispersive in space, it is not very likely to compose a tight object. To compensate the spatial 
dispersion, one can introduce some form factors at the effective vertices Q. 

On contraries, we can be convinced that if there are two heavy quarks {bb, be, cc) in the 
baryons, they would tend to constitute a substantial diquark with small spatial dispersion which 
can serve as a static 3 color source for the light quark P]P|. Savage and Wise estimated spectra of 
baryons with two heavy quarks in the heavy quark effective theory (HQET) R. Recently, Ebert 
et al. evaluated the spectra of such baryons in terms of the local Schrodinger-like quasipotential 
equation [^. In the framework of the potential model, the interaction between the light quark 
and heavy diquark can be derived by calculating their elastic scattering amplitudes Q, but 
the key point is the form of the effective vertices at the diquark-gluon interaction. Similarly, 
Gershtein et al. also considered the spectroscopy of doubly charmed baryons where they include 



angular and radial excited states [10| 



As aforementioned, even though the diquark consisting of two heavy quarks is tight, in 
reality it is still not point-like and the small deviation should be estimated. The authors of 
ref.Q introduced simplified form factors and ignored their /c^— dependence to take the effects 
into account. 

In this work, we re-derive the effective potential by using the B-S equation and obtain the 



effective vertices. In the derivations, the fc^— dependence is retained exphcitly. We find that this 
dependence leads to an extra Yukawa-type term. 

Because of the serious relativistic effects of the hght quark, the non-relativistic expansion of 
the potential becomes dubious, since we must truncate the expansion, usually to order p ^/m^, 
where m^ may be the light quark mass. Generally, p^ ~ ^qcd ~ (^-^ GeV)'^ , m takes the 
constituent quark mass as m„ ~ rriii ~ 0.33 GeV, thus the expansion factor v ^/c^ = p "^ /m? ~ 
0.37, so the next leading order should be of certain contributions. In fact, one can attribute such 
uncertainties to the parameters which exist in the potential, but cannot be directly measured. 
Therefore considering this situation, instead of the T and J/^ spectra we re-fit the data of 
S'*)— meson, which includes a light quark and a heavy b-quark, to obtain the parameter values. 
Substituting the re-obtained parameters into the potential we evaluate the baryon spectra, one 
has reason to expect that in this way, the errors can be substantially reduced. 

In the process, we also consider the operator-ordering problem. Since we derive the scattering 
amplitude in the momentum space, all quantities are commutative. But when the Fourier 
transformation with respect to the exchanged momentum k is carried out, the coordinate r and 
momentum p co-exist, their ordering becomes a problem. This is an inconsistency which exists 
in the treatment. However, our results (see below) indicate that the ordering only determines 
the parameter values, but the final measurable quantities, i.e. the spectra do not deviate much 
from each other in different ordering schemes. 

We also briefly discuss the possible mixing problem and point out the drawbacks of the 
potential model. 

The paper is organized as following. After this introduction, we present the formulation and 
explicit form of the potential. We also discuss some concerned problems such as the ordering 
schemes. In Sec. Ill, we discuss how to obtain the spectra in terms of variational method, trial 
function and concerned issues. In Sec. IV, we give the numerical results and the adopted param- 
eters. The last section is devoted to our conclusion and discussion. 



II. Formulation 

(a) The effective vertices for diquark-gluon coupling. 

Since the diquark is not rigorously a point-like subject, we cannot simply use the vertices 
in the fundamental QCD theory. Instead, we derive such effective vertices with certain form 
factors in terms of the B-S equation. The diquark contains two heavy quarks which constitute 
a color 3 triplet, for this bound state the Cornell potential would be a good approximation and 
we use it as the B-S kernel. 

We derive the effective vertices for SSg, AAg and SAg, ASg as following: 

{S'iv')\J,,\Siv)) = VmAPUW^ + f2V^) for SSg coupling, (1) 

+ MV v'){v'* ■ v)v'^ + Mri- v'){r,'* ■ v)v^ 
+ frv'/liv v') + fsiv'* ■ v)ri^ 

+ fioief,upaV'*'"v^v''] for AAg coupling, (2) 

{A'{v',r,')\J,\S{v)) = VMM^[hW; + fi2{v'*-v)v'^ 

+ fi^iiv'* ■v)v^ + fiiie^yp„r}'*"v'Pv''] for ASg coupling, (3) 

where S and A stand for scalar and axial-vector diquarks, v' ,v,rj' ,r],M' ,M are the four- 
velocities, polarization vectors (for axial- vector diqaurk only), and masses of the diqaurks in 
the " final" and " initail" states of the scattering respectively. The corresponding form factors are 
derived by solving the transition B-S equation and the details were given in our previous work 
1^]. In our case, we find relations 

h = f2 = h = h = -h = -h = /i4 = /, (4) 



and 

/g = /lo = /ii = /i2 = /i3 = 0, (5) 

and the relations (|5|) can be realized by simple parity analysis. Obviously, the terms related 
to /s and /g are proportional to |vp (Ipp/m-^) so that can be neglected as we only keep the 
non-relativistic expansion up to order of p ^ . 

Here we would like to draw attention of readers that in expressions (pi) and (pi), the order 
of 7] and rj'* is not trivial. When we derive these formulae in Quantum Field Theory (QFT), 
they are commutative, so that we can put them in any order, however, as we turn r] and rj'* into 
spin-operators of Quantum Mechanics (QM), the order problem emerges. Because S— operator 
is not self-commutative, thus in the operator form, r] ■ r]'* ^ ry'* • rj. Therefore, when we write 
the expressions, we must be very careful about the order. Our strategy is that we keep the right 
forms of the leading terms relating to spins such as the spin-orbit coupling terms to choose the 
appropriate orders. In fact, we can also determine the order by analyzing the symmetry of the 
whole Hamiltonian as we have done in the expressions (^) and (^). 

The form factors /» involves the B-S integrals and cannot be analytically expressed. One can 
only obtain the numerical results instead. However, in order to serve our final goal to derive an 
effective potential, we need an analytical expression for the Fourier transformation. So we have 
simulated the numerical results with various function forms, finally we decide 

f _ A B 

k2 ~ k2 + k 2 + C2 ' ^^' 

where k is the exchanged three-momentum, gives the best fit. The parameters A,B and C are 
numerical values. In this expression, we keep the explicit A;— dependence of the form factors. In 
fact, the expression (^) can be rewritten as 

_l_ ^ iA + B)ik^ + ^C') 

k2 k2(k2 + C2) ' ^^ 

and k2 -|- 6*2 = — (A;2 — C*^) at the case ko = 0. It is the familiar pole-like form factor which is 
widely used in phenomenology Ef. 



(b) We derive the effective potential by calculating the elastic scattering amplitude, then 
we need to turn the corresponding quantities into the quantum mechanics operators. The 
polarization vectors 7] or rj' of the axial vector diquarks must be normalized as r/^ = r/ ^ = — 1 
according to the quantum field theory. Turning ??(r/') into QM spin-operator, we have 

V = ^{0-Sh,S + ^{P-S)P), (8) 

where S is the spin-operator, (3 = p/M, 7 = E/M are the boost factor and p, E^ M are the 
momentum, energy and mass respectively, the factor l/\/2 guarantees the right normalization 
for the axial- vector diquark S{S + 1) = 2. 

In derivation of the potential, we first calculate the part corresponding to the scattering 
amplitude induced by one-gluon exchange, 

M5^«°"(p,k) =< A'^A'^ > gl ^ u(p;)7'^n(pi)I)^.(k)(p^|J^|p2), (9) 

lQ^EiE[E2E'^ 

where the Coulomb gauge for the gluon propagator is chosen, 

1 1 h^y 

^°°(^) = -^^D^'ik) = -^{5'^ - -^),D'^' = i^^o = 0. (10) 

Thus we have the expressions for the transition amplitudes as following: 

..qq / , N IGiras f r p^ k^ , 1 1 , 

F,f„„(p,k) = ^^[1 + -^^ ( +_) 

^ 3 k^ 7711771,2 4r77i Znii 771,2 

iSi • (k X p) 1 2 

+ ^ ( + )], (forSq^Sq); (11) 

zr7ii nil 1^2 

^ 3 k^ 77711712 4r7ii 2r7ii 1712 

^ iSi-(kxp) 1 ^ 2 iS2-(kxp) 1 ^ 1 
2i7ii mi 1112 Am2 nii m2 

(S..S,)k=-(S,-k)(S^|^ (faAq^Aq); (12) 



4l71ir7i2 

16 7ras / 7S2(k x p) 1 J_ 
3 2\/2k2 m2 mi 1712 

Si-S2)k2-(Si-k)(S2-k) 



ni"o„(p,k) = ^"r.io [^^^^^^^(z^ + 



77117712 



], (for Aq ^ Sq or Sq ^ Aq). (13) 



The Fourier transformation would bring up the ordering problem which will be discussed below, 
and then we will present the expressions in the configuration space in subsection (e), while the 
confinement part is given in subsection (d). 

(c) Ordering of operators. 

When we derive the scattering amplitude in the momentum-space, all quantities are com- 
mutative, however, when we transform them into the QM operators and carry out a Fourier 
transformation to the configuration space, there exists an ordering problem in general. 

When we transform k(= P]^ — Pi) into r, which in fact is the relative radial vector between 
the light quark and the heavy diquark, if we choose the center of mass of the system to be the 
coordinate origin, r2 of the heavy diquark would be very close to zero and ri r^ r. The reduced 
mass of the system mM/{m + M) -^ m, which is almost the mass m of the light quark. The 
momentum p would remain as a derivative operator in the configuration space, thus the ordering 
problem emerges. 

For example, there are three different orders for p, p,g{r) where g{r) is a function of r, as 

5(r)p^ p-g{r)p, -[g{r)p'^ + 2p ■ g{r)p + p'^g{r)]. 

In most literatures, one just simply takes g{r)p ^. In our work, we compare the different 
ordering schemes and our numerical results show that different schemes would lead to different 
parametrizations, but the final measurable spectra are not sensitive to the ordering at all. 

(d) The confinement part. 

The confinement part of the potential is fully due to the non-perturbative QCD effects and 
is not derivable in any established theoretical framework. So far, one can only postulate its form 
and determine the concerned parameters by fitting data. 

The most commonly adopted confinement form is the linear potential Vconf = ar + b at 
the leading order. It can be split into a scalar and vector pieces which may lead to different 
relativistic corrections. Since its source is obscure so far, one cannot decide fractions of each 



piece. But in general, it can be written as |g] 

14on/(0 = ^(ar + 6), (14) 

con J \ 



VclUr) = {l-K){ar + b). (15) 



The resultant potentials with all relativistic corrections are 

VLfir) = yco'nfir)-^P-V'%r)p + ^^AV^%r) 



1 V'^ir) 1 V'^jr) ^ 

-Si -L- 7^-2 Z Ss-L, (16) 



2mf r 2m2 r 

•' ■' Ami 2/U 2m2 mi 



m,im,2 m-i fi Zmi r 



1 ^n^ s^ j^ , 2(1 + K 

2m2 r 377117712 



S2 • L + ^ Ay^"(r)Si • S2 



+ :^i^^^^^-V"^\r))Si„ (17) 

077717772 T 



and 



3 

512 = Si • S2 2 (^1 ■ ^)i^'2 ' ^)- 

In later numerical calculations, we choose several values of k. 

(e) The potential. 

Finally we have the full Hamiltonian 

H = K + V, (18) 

where K is the kinetic part and 

V = Vgluon + Vconf. (19) 

and 

yconS = Vconf + V^conf- (20) 



The single gluon exchanged potential Vgiuon has the following forms as 

T/'S'S' / ^ _ 
^ gluonV J 



4a., A 



tAA 



ySA 



gluonx ) 



„ . , , -p2 + 2p--p + p'- 



mi 2mi 7712 2mi mi m2 ?''^ 



+ 



p-Cr 



-Cr 



-Cr 



+ 



, r + 2p p + r 

r 4?7iim2 r r 

BC\ 1 1 .e-'^^ 

+ 



4mi 2mi m2 r 

B , I 2 ,(Cr + l)e-^'' 

( 1 ) 

2?7ii mi m,2 r^ 



Si • L], (for scalar — diquark + q baryons);(21) 



^ gluonV I q l ' 



A 1 



1. -21. 



+ 



+ 



q L . . ,-p" + 2p--p + P^- 

6 r 4m,im,2 r j. j. 



A , \ 2 ,Si-L 

(— + 



mi 2mi m,2 2m,i mi ra2' r 

A Si2 27rA 

-. 3" - ^ (Si • S2)5(r) 

4mim2 r-' om,im,2 



^.1^1 . S2-L 
4m2 mi m2 r^ 



Be 



-Cr 



+ 



B 



,-Cr 



,-Cr 



-Cr 



r A77li77l2 r 

BC^ , 1 1 .e-'^^ 



( — r + 2p 



-p + p 



-(TTI- + — 



4mi 2m,i m2 t 2m,i m-i m2 
H ) 5 ^2 • -L 



5,12 (Cr + l)e-^ 

1 j ^ Si • L, 



4m2 mi m,2 r"' 

B C^r^ + 3Cr + 3 



12m,im,2 



e-^"5i2 



B C'e 



2^-Cr 



6mim2 r 
4ac, A 



-Si • S2], (for axial — vector — diquark + q baryons); (22) 



1 1 ,S2-L 



A S 



12 



+ 



AttA 



3 2\f2ni2 "rn\ 7712 r^ 2\/2mim2 r^ 3\/2mi 

B I I (CV+2)e^„ ^ 

(, \ j n =>2 • L, 



m2 



(Si • S2)<5(r) 



2\/2m2 mi 7712' 

B CV + 3Cr + 3 _cr 



6v2W'l'722 T 

B C^e-^'- 



5 



12 



Si • S2], (for mixing between scalar — diquark + q 



3\/2mim2 r 

and axial — vector — diquark + q baryons), 



(23) 



where 

Su = ^^"^-y^-^^ - S, . S. (24) 

In the expressions A and S stand for the axial-vector and scalar respectively. Later we will 
show that even though Vf^„ derived in QFT is not trivially zero, in the non-relativistic QM 
framework, it can give only null contribution. We will discuss this issue in the last section. 

III. The variational method 

(a) We choose the variational trial function with a single parameter for the IS state as 

i;ir,9,^)=R{r)Yoo{9,^), (25) 

and 

R{r) = Ne-^''\ (26) 

where 5 = 4/3 and the normalization is 






and A is the variational parameter. This form is discussed in our earlier work |11] where we 
tested some 5— values and found that 6 = 4/3 would be more appropriate for the Cornell-type 
potential, we will discuss this problem further in the last section. 

(b) The parameters as, a and b. 

As noticed, the relativistic effects are serious because of the existence of a light quark. Unlike 
the heavy quarkonium, such as J/ip,T etc., truncation of the non-relativistic expansion where 
we only keep it up to p "^ /m?' order, is not a good approximation. However, we can partly 
compensate the effects by attributing the uncertainties to the potential parameters which are 
not directly measurable. In other words, in the process of fitting data of mesons containing a 
light quark, such as B^*\ we have attributed the unknown factors into the phenomenological 

10 



parameters, then later when we use the set of parameters to evaluate the spectra of baryons 
containing two heavy quarks and a light quark, the non-perturbative QCD effects and the 
relativistic influence are or at least mostly included. Obviously in the case, if one used the 
parameters obtained by fitting data of heavy quarkonia, the errors are un-controllable. Here we 
choose B^*' data to obtain as,a and 6. It is worth noticing that in the D-case, the relativistic 
effects are serious and the charm-quark is not heavy enough, so when we apply our trial function 
to the D-case, we find the minimum of the expectation value of energy is not stable. Thus we 
abandon the D-case. When we use the variational method to obtain the parameters, we retain 
all the relativistic corrections in the potential for i?^*^— mesons. 

(c) Then we turn to calculate the spectra of the baryons containing two heavy quarks, thus 
A stands as the variational parameter. The expectation value of H is 

where -R(A) is the chosen trial function (p^). Then minimizing E{X) as 

dE{Xl_ 
dX 

we obtain the A— value. In the expression H is the full Hamiltonian given in eq.([l8|). 

The advantage of using the variational method is obvious, that is we are able of treating all 
terms simultaneously. Unlike the perturbation method where all relativistic corrections which 
are very large in this case are dealt with perturbatively, so that remarkable errors for the baryons 
which contains not only two heavy quarks but also a light one, emerge due to the ill-treatments, 
by contraries, the ambiguities can be avoided in our treatments. 

IV. The numerical results 

We are listing some concerned parameters which appear in our formulae. 
The A, B, C— values in eq.(^) are 

A =1.0060428, 5 = -1.280532, C = 2.425356, for cc diquark; 

11 



A = 1.0019211, B = -1.545131, C = 6.415003, for be diquark; 
A = 1.0005426, 5 = -1.244000, C = 3.676031, for bb diquark. 
We have the constituent quark masses and the heavy diquark masses as 



m„ 



md = 0.33 GeV, nis = 0.5 GeV, Mcc = 3.26 GeV, Mbc = 6.52 GeV, Mbb = 9.79 GeV. 



It is noted that bb and cc diquark must be axial vectors, but be can be either a scalar or an axial 
vector, the mass splitting of the scalar and axial- vector be diquarks can be neglected in practical 
calculations. 

The baryon spectra are calculated and the results are given in Tables 1 and 2, with q being 
u or d. In Table 1, we choose k = —1 for the confinement potential (^) which is consistent with 
that used in ref.|^], and list results corresponding to various ordering schemes. In Table 2, we 
change the k— values in the confinement potential and use the ordering scheme 2, i.e. p • g{r)p. 

Table 1. 



type 


Ordering 1 


Ordering 2 


Ordering 3 


(ccq)(l/2) 


e=0.1972, Mb=3.787 


e=0.2175, Mb=3.808 


e=0.1480, Mij=3.738 


(ccq)(3/2) 


e=0.2229, Mb=3.813 


e=0.2425, Mb=3.832 


e=0.1793, Mij=3.772 


(ccs)(l/2) 


e=0.2027, Mb=3.963 


e=0.1101, Mb=3.870 


e=0.0875, Mb=3.851 


(ccs)(3/2) 


e=0.2254, Mb=3.985 


e=0.1375, Mb=3.897 


e=0.1187, Mb=3.879 


(cbq)(l/2)5 


e=0.2411, Mb=7.091 


e=0.2325, Mb=7.082 


e=0.1931, Mij=7.043 


(cbq)(l/2)A 


e=0.2320, Mb=7.082 


e=0.2234, Mb=7.073 


e=0.1823, Mb=7.032 


(cbq)(3/2)A 


e=0.2455, Mb=7.095 


e=0.2369, Mb=7.087 


e=0.1981, Mij=7.048 


(cbs)(l/2)5 


e=0.2268, Mb=7.247 


e=0.1202, Mb=7.140 


e=0.1158, Mij=7.236 


(cbs)(l/2)A 


e=0.2181, Mb=7.238 


e=0.1096, Mb=7.130 


e=0.1073, Mij=7.225 


(cbs)(3/2)A 


e=0.2310, Mb=7.251 


e=0.1253, Mb=7.145 


e=0.1210, Mij=7.241 


(bbq)(l/2) 


e=0.2168, Mb=10.337 


e=0.1983, Mb=10.318 


e=0.1731, Mb=10.293 


(bbq)(3/2) 


e=0.2251, Mb=10.345 


e=0.2069, Mb=10.327 


e=0.1830, Mb=10.303 


(bbs)(l/2) 


e=0.2096, Mb=10.500 


e=0.0910, Mb=10.381 


e=0.1061, Mb=10.396 


(bbs)(3/2) 


e=0.2174, Mb=10.507 


e=0.1009, Mb=10.390 


e=0.1150, Mb=10.405 



In Table 1, "ordering 1" means g{r)p'^, where Og = 0.23, a = 0.11,6 = —0.13; "ordering 2" 
means p • g{r)p with Os = 0.41, a = 0.09, b = —0.21; "ordering 3" means {g{r)'p ^ + 2p • g(r)p + 
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p ^(7(r))/4 with as = 0.23, a = 0.11, h = —0.27. The subscript A and S stand for the axial vector 
and scalar respectively, e is the binding energy and Mb is the baryon mass with unit GeV. In 
the calculations, k = — 1. 



Table 2 



Type 


K = 


K = 0.5 


K= 1.0 


(ccq)(l/2) 


e=0.2133, Mb=3.703 


e=0.2265, Mb=3.817 


e=0.1713, Mij=3.761 


(ccq)(3/2) 


e=0.2449, Mb=3.735 


e=0.2583, Mb=3.848 


e=0.2153, Mb=3.805 


(ccs)(l/2) 


e=0.1466, Mb=3.807 


e=0.0834, Mb=3.843 


e=0.0501, Mb=3.810 


(ccs)(3/2) 


e=0.1756, Mb=3.836 


e=0.1180, Mb=3.878 


e=0.0960, Mb=3.856 


(cbq)(l/2)5 


e=0.2539, Mb=7.10A 


e=0.2582, Mb=7.108 


e=0.2294, Mij=7.079 


(cbq)(l/2)^ 


e=0.2428, Mb=7.093 


e=0.2463, Mb=7.096 


e=0.2142, Mij=7.064 


(cbq)(3/2)^ 


e=0.2593, Mb=7.109 


e=0.2640, Mb=7.114 


e=0.2368, Mb=7.087 


(cbs)(l/2)5 


e=0.1799, Mb=7.200 


e=0.1045, Mb=7.125 


e=0.1048, Mij=7.125 


(cbs)(l/2)A 


e=0.1693, Mb=7.189 


e=0.0897, Mb=7.110 


e=0.0885, Mb=7.108 


(cbs)(3/2)A 


e=0.1851, Mb=7.205 


e=0.1116, Mb=7.131 


e=0.1126, Mb=7.133 


(bbq)(l/2) 


e=0.2131, Mb=10.333 


e=0.2187, Mb=10.339 


e=0.1748, Mb=10.295 


(bbq)(3/2) 


e=0.2236, Mb=10.344 


e=0.2298, Mb=10.350 


e=0.1894, Mb=10.309 


(bbs)(l/2) 


e=0.0876, Mb=10.378 


e=0.0677, Mb=10.357 


e=0.0607, Mb=10.351 


(bbs)(3/2) 


e=0.0997, Mb=10.390 


e=0.0813, Mb=10.371 


e=0.0759, Mb=10.366 



In Table 2, for the confinement potential Vf^^j = n{ar + b),V^^f = (1 — K){ar + 6), we 
have: for k = 0, the fitted ag = 0.44, a = 0.14,6 = —0.37; for k = 0.5, the fitted as = 0.5, a = 
0.14,6 = —0.37; for k = 1. the fitted Us = 0.73, a = 0.16,6 = 0.45. Here we use the ordering 2, 
i.e. p-g{r)p. 



V. Conclusion and discussion 

As we noted in the introduction, we calculate the spectra of baryons containing two heavy 
quarks which can constitute a diquark. It is believed that such a subject can be spatially tight 
and serve as a color source for the light quark. This picture greatly simplifies the calculations. 
The difference of our method from previous works is in several aspects. 

First, we use effective vertices DD'g where D and D' are scalar or axial- vector diquarks and 
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g is gluon. We derive the form factors at the vertices based on the B-S equation and so we can 
keep their exphcit A;^— dependence which leads to an extra Yukawa- type term in the potential. 

Secondly, we investigate the ordering problem which is brought up by the Fourier trans- 
formation with respect to the exchanged momentum k and the quantization of momentum p. 
We find that various ordering schemes can lead to different parametrizations for as (note that 
here a^ is an effective coupling constant, but not that from the fundamental QCD theory), and 
a, b which are not directly measurable. However, we find that the final results do not deviate 
much from each other. So we can conclude that at least the ordering schemes do not seriously 
influence the spectra evaluation in the variational method. In our future work, we will continue 
to investigate if the ordering schemes can induce other observable effects such as the effective 
decay constant etc. 

In our scenario, we use the variational method and the Hamiltonian includes not only the 
leading kinetic and potential terms, but also the relativistic corrections up to order of p^. 
Because the relativistic effects are very serious in the case where a light flavor is involved, this 
treatment is superior to the perturbative method. To reduce the uncertainties and errors brought 
up by the truncation of the non-relativistic expansion, we use the B^*' data where a light quark 
is moving around the heavy b-quark, as inputs to obtain suitable parametrization. As hoped, 
most of those uncertainties and errors can be attributed into the phenomenological parameters 
Us, a and b. 

As well-known, the potential model cannot perfectly describe hadron characteristics which 
are mostly determined by non-perturbative QCD effects and we have no reliable knowledge on 
it so far. But as long as we use the experimental data as inputs to parametrize the model, the 
disadvantages can be partly compensated. In the trial function i?(A, r) in the form of eq.(p5|), 
we priori take 5 = 4/3 based on our previous work. In fact, 5 should be an irrational value 
between 1 and 2 which correspond to the solutions for the Coulomb and harmonic oscillator 
potentials. But 6 = 4/3, as indicated in ref.[^, is a satisfactory value for the Cornell potential. 
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Numerically, in our previous work[ll|, we obtained the irrational number 6 for the Cornell 



potential as 5 = 1.33809..., which is very close to 4/3. Then we take the value and only let A 
be the unique variational parameter, this treatment can greatly simplify the calculations with 
sufficient accuracy being kept. 

Even though the diquark picture is believed to work in this case and the derived form 
factors further improves the situation, there still exists small deviation from reality, including 
the diquark masses. This should be further investigated. 

Finally, as we pointed above, although the mixing term in eq.(|23|) which is derived in QFT 
is not trivially zero, when we sandwich it among the quantum states, we have 

< V(l/2, A, I = 0)\V^\l'^lmi/2, S,l = 0) >=< ^(1/2, A, I = l)\V^\ltlmi/2, S,l = l) >^ 0. 

(28) 
The matrix elements are absolutely zero. The reason is simple, because in the framework of non- 
relativistic quantum mechanics, there are no creation and annihilation operators as in QFT, so 
that we can only deal with elastic scattering. The mixing between -0(1/2, A) and ^^(1/2, S) refer 
to a change of spin or particle identity, so cannot appear in QM even though we know such mixing 
must exist and may play important roles to hadron spectra. For example as in a completely 
different area of the hadron spectroscopy, the mixing between glueball and quarkonium is known 
as very important or even crucial to phenomenology, but we cannot evaluate it in the potential 
model. We will further study these mixing effects in our future work ^ 



The B-factory and other facilities of high energy experiments may provide data on Sec and 
other such baryons. Once the data are available, we may re- adjust our input parameters and 
make further predictions on the spectra and other characters of the baryons, then we can testify 
the validity of the diquark picture and the non-relativistic potential model. 
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